function FC = FWC(x)
    global rho rhopl K L;
    c=x; b=1/sqrt(1-rho); a = (rhopl-(b+L*c)*rho)/(1+(K-1)*rho);
    ret=K*a*(a*(1+(K-1)*rho)+2*(b+L*c)*rho)+(2*b+L*c)*c*(1-rho)+rho*(b+L*c)^2;
    
FC = ret;